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We study oscillons, extremely long-lived localized oscillations of a scalar field, with three different 
potentials: quartic , sine-Gordon model and in a new class of convex potentials. We use an absorbing 
boundary at the end of the lattice to remove emitted radiation. The energy and the frequency of an 
oscillon evolve in time and are well fitted by a constant component and a decaying, radiative part 
obeying a power law as a function time. The power spectra of the emitted radiation show several 
distinct frequency peaks where oscillons release energy. In two dimensions, and with suitable initial 
conditions, oscillons do not decay within the range of the simulations, which in quartic theory 
reach 10 8 time units. While it is known that oscillons in three-dimensional quartic theory and sine- 
Gordon model decay relatively quickly, we observe a surprising persistence of the oscillons in the 
convex potential with no sign of demise up to 10 7 time units. This leads us to speculate that an 
oscillon in such a potential could actually live infinitely long both in two and three dimensions. 
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I. INTRODUCTION 

The theorem by Derrick states that there are no non- 
trivial static solutions of finite energy in scalar theories 
above one dimension [l|. This does not say anything 
about time-dependent solutions. Example of such are Cu- 
bans , whose stability is guaranteed by a conservation 
law (for more recent studies on Q-balls see e.g. [HIH) and 
that have been studied e.g. as a dark matter candidate 
(see for instance Q). There is, however, no obvious rea- 
son why oscillations of a real scalar field can persist and 
stay localised, yet such configurations, oscillons exist in 
many models (we will use the term oscillon, presented 
originally in [6j, in a slightly loose manner throughout 
this study). There is a stable breather solution in one- 
dimcnsional sine-Gordon model that during its period de- 
forms to a separated kink and anti-kink. There is no such 
exactly stable solution in 4 theory 0. However, even 
though the solution radiates energy away, there remains a 
possibility that the radiation rate is so suppressed that it 
will take indefinitely long for the dissipation to complete. 

Oscillons were first found already in the 70's 0, 
and then rediscovered in the 90's when the dynamics of 
phase transitions was studied [||. There has been sub- 
stantial interest in oscillons recently. On the theoretical 
side of understanding the dynamics of oscillons, there 
have been a number of studies of oscillons of 4 theory 
in small-amplitude approximation [Tol - [l3j . Furthermore, 
the attractor basin of oscillons and its fractal nature have 
been studied both analytically [l4| and numerically [l5j |. 
Very recently there has been interest in oscillons coupled 
to gravity, oscillatons ^6l - [l~8j . Oscillons have also been 
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found in dilaton-scalar theories [19j. A new class of so- 
lutions, called flat-top oscillons by the authors, were re- 
ported in [2^ . While the body of work mentioned above 
has been on a classical level, the study in [2l[ considered 
quantum corrections oscillons are subject to, redeeming 
oscillons considerably less robust against the quantum 
effects than often previously has been assumed. 

Though stability or finite lifetime are of considerable 
theoretical interest, from the point of view of a many 
phenomcnological consequences, it is not important if the 
non-trivial solutions are actually stable like the breather 
in one-dimensional sine-Gordon model as long as they 
are far longer-lived than the natural time scale related to 
the process on which they will have an effect. A natu- 
ral realm of oscillons to appear is in the early Universe 
and such a process there could be baryogencsis. The role 
oscillons could potentially have on baryogenesis is provid- 
ing the necessary non-equilibrium needed for successful 
creation of matter-antimatter asymmetry analogous to 
Q-balls whose cosmological impact has been studied ex- 
tensively (for a review, see e.g. [22|). In this context it is 
very interesting that an oscillon solution has been found 
in the bosonic sector of the Standard Model (23rl25j ). 
When the scalar and vector masses in the theory are set 
to be inn — 2mw, the oscillon is sufficiently long-lived 
that it has not been seen to decay in any numerical simu- 
lation to date (23rl25| . Oscillons have also been found in 
two and three-dimensional Abclian-Higgs model in deep 
type I regime [2(1 [27j . A number of studies have consid- 
ered oscillons in the early Universe. While the oscillons 
in the Standard Model or Abclian-Higgs model are obvi- 
ously objects in several fields, a recent study [28[ found 
that oscillons can also form in two scalar fields in mod- 
els that are relevant from the point of view of hybrid 
inflation (on that topic see also (29j|)- Other recent stud- 
ies of formation of oscillons after an inflationary epoch 
but in single field models include [10, HH . Regarding the 
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subsequent evolution of oscillons after their formation in 
a cosmological setting, both numerical [32[ and analyti- 
cal [33[ considerations in one dimension have shown that 
oscillons can persist for a considerably time in expanding 
backgrounds (see also Hq|). The proposal of oscillons fa- 
cilitating vacuum tunneling in the context of the string 
landscape [34| is an application quite far from the original 
idea of oscillons affecting the bubble nucleation process, 
the resonant nucleation studied e.g. in [35J (for a study 
of bubble collisions, see also |3{|). A recent study consid- 
ered the formation of oscillons in the quintessence field 
at the present epoch of the Universe |37| . 

Lastly, regarding the significance of the oscillons, it 
should be mentioned that they are an interesting link 
between particles and solitons. Several numerical stud- 
ies in one dimension have considered the formation of 
kink-antikink pairs. It has been found that scattering of 
breathers can create a kink and an antikink [38| . Os- 
cillons can also decay to a kink and antikink when the 
potential of the model is distorted [39] , as well as oscillons 
being an intermediate state in the process of formation 
of the kink-antikink pair |4(j ■ 

In an earlier work (4l| we studied oscillons in two di- 
mensions on the lattice with periodic boundary condi- 
tions. The aim of this study is to examine in detail the 
evolution of oscillons and the properties of the radiation 
they emit, in order to better understand the reasons for 
their longevity with the use of absorbing boundary con- 
ditions. 

Previous work in a theory with a quartic potential has 
identified stable periodic solutions with incoming radi- 
ation, called quasl-breathers, which are closely related 
to oscillons [42|, |43|. There is a critical oscillation fre- 
quency at which the quasi-breather has the minimum en- 
ergy [42l. [43l] , and once the incoming radiation is removed, 
the frequency of the resulting solution evolves towards 
its critical value according to a power law in time [43j . 
Power laws featured also strongly in [3, Ell where the 
radiation rates of oscillons were studied, starting from 
the assumption of a strictly Gaussian form for the oscil- 
lon, and a decay width calculated by comparison to linear 
theory. An equation relating the energy loss rate to the 
rate of change of the oscillon amplitude was derived, with 
power law solutions for the time evolution. Note that os- 
cillons in the small-amplitude expansion [l2|, [l3| behave 
differently, in that they shed energy in an exponentially 
suppressed way as the amplitude goes to zero, and do not 
obey a simple power law. 

In this paper we determine these power laws in three 
models in two and three dimensions, consisting of a sin- 
gle canonically normalised real scalar field <p with poten- 
tials listed in Table HI and examine the power spectra 
of the emitted radiation. Our power spectra show that 
for the quartic and sine- Gordon theories the radiation is 
predominantly emitted at an integer multiple of the ba- 
sic oscillation frequency, either 3 or 2 depending whether 
the potential is symmetric or not. For the convex po- 
tential, the radiation is emitted just above the threshold 



model 


potential V((j>) 


quartic (0 4 ) 
sine- Gordon 
convex 


— (1 + C0S(7T(^>)) 
T 

1 + (ft 2 ? 



TABLE I: Potentials for the three models considered in this 
paper. In the convex potential, p is a real parameter, with 
p < 1. For the data shown p = 0.45. 



set by the mass parameter. A fundamental assumption 
made in (45j in developing a theory of oscillons in the 
quartic theory was that the radiation is predominantly 
emitted just above threshold, which we see is incorrect. 
Nonetheless, good fits to the time evolution of the energy, 
frequency and amplitude were obtained in that work. 

An important technical development in our work is 
the construction of absorbing boundary conditions for 
a massive field, which provides an economical alternative 
to the adiabatic damping technique, introduced in [46| . 
and also used in (43|, which sets a location dependent 
friction co-efficient on the lattice. Absorbing boundaries 
provide an alternative approach to remove the dispersive 
waves from the lattice, requiring no extra lattice sites for 
its operation. 

The paper is organised as follows. The numerical set- 
up is briefly reviewed in the following section. While 
the equations of the absorbing boundary conditions are 
presented in the appendix we briefly discuss their impli- 
cations for simulations together with the data obtained 
using quartic theory as an example. We then proceed to 
present the results for the power laws and the radiation 
power spectra in the three theories in two and three di- 
mensions. We find a surprising stability for the oscillon 
in the convex potential, which leads us to conjecture that 
it may be stable in three dimensions. 



II. NUMERICAL SET-UP 



by 



The Lagrangian for a single real scalar field 4> is given 



(i) 



and the equation of motion thus reads 

4>- V 2 <P + V'{(f>) = 0. (2) 

This field equation is evolved using a leapfrog update 
and a three-point spatial Laplacian accurate to 0(dx 2 ). 
We perform (l+l)-dimcnsional simulations assuming ra- 
dial symmetry These permit the choice of a much finer 
lattice spacing compared with (2+l)-dimensional simu- 
lations carried out in a previous study [4l| . We kept the 
ratio of time step and lattice spacing fixed and match- 
ing to the one used in (2+l)-dimensional simulations, 
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dt : dx = 1 : 5. The size of the onc-dimcnsional grid was 
chosen so that the physical distance from the core of the 
oscillon to the absorbing boundary was always 100 units. 
The simulations are performed on a grid of 10,001 lattice 
points and dx = 0.01, dt = 0.002. 

We set the field <j> an d the field momentum II = <j> at 
the origin r = to have the same value as on the next 
lattice site with a non-zero value of the radius r. The up- 
dates used at the far end of the grid, where the absorbing 
boundary is located, are derived in the appendix, while 
we comment the simulations with the absorbing bound- 
aries shortly in the next section. We use three consecu- 
tive points in the grid to evaluate the gradients with the 
precision also accurate to 0(dx 2 ). 

III. QUARTIC POTENTIAL 

We study here the degenerate double-well quartic po- 
tential which defines the model that we refer also as 4 
theory 

y(0) = iA(^ 2 -« 2 ) 2 , (3) 

which sets the mass of the theory to be m 2 = 2Xv 2 . 
Scaling out the vacuum expectation value and coupling 
the potential can be written 

V{^)^\{^-lf (4) 

so that the minima are located at cf> = ±1 and the local 
maximum at <p = 0. We use the following Gaussian initial 
data to create an oscillon 

0(r)= (l-C-cxp(-r 2 /r 2 )), (5) 

where r is the distance to the center of an oscillon. The 
width of the distribution is set to be r ~ 2.9 (in units 
of (y/Xr))" 1 ) in two and tq ~ 3.0 in three dimensions, 
suggested optimal choices in j47|. The maximum dis- 
placement in two dimensions is set to be C = 1 so that 
the center of the oscillon starts from the local maximum 
of the potential, while we use C = — 1 in three dimensions 
again following [47| ■ 

A. Properties in (fi 4 theory in two dimensions 

As we kick off the oscillon by the somewhat arbitrary 
initial condition ([5]), it naturally sheds at least some 
fraction of the energy in the form of propagating waves 
moving away. As long as that emission is not too vi- 
olent and at least further distance away from the os- 
cillon core spreading and consequently damped in two 
or higher dimensions these waves can be assumed to be 
well-characterized just by a linear approximation of the 
equation of motion ([2]) , see (|Al[) . Once the wave reaches 
the boundary of the grid it is absorbed and the energy 
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FIG. 1: The total energy in the lattice (top) and energy in- 
side a shell of radius r = 5ro (bottom) in <j} 4 theory in two 
dimensions. 



carried by the wave removed from the lattice. We thus 
expect monotonically decreasing total energy in the sys- 
tem, and this is also what we observe. Figure [1] shows 
the total energy (top line) in the lattice at the very early 
stage of the simulation. The total energy in the lattice 
will always include the radiative component emitted by 
the oscillon which has not yet reached the boundary of 
the grid. Thus we measure also the energy inside the ra- 
dius r = 5ro from the center of the oscillon and refer this 
as the energy of the oscillon. The oscillon is generally well 
localised inside this volume, but it covers less than 15% 
of the lattice used here. The energy inside this shell is 
also depicted in Figure [T] (bottom line). In the beginning 
it takes a finite time the wave to reach the boundary and 
even then the decrease in total energy is more gradual 
compared with the energy inside the shell; this is caused 
by the dispersion of the emitted waves. However, after 
the initial transient phase these two ways to measure en- 
ergy track each other well; there is only a time delay 
between them. This also provides a crucial check for the 
adequacy of the absorbing boundary conditions we use. 
Any reflected waves will travel back and reach the inner 
shell consequently increasing energy measured inside the 
radius there. There is no sign of such a burst at a vis- 
ible level in Figure [T] and we expect only a very minor 
reflection taking place at the boundary. We monitor the 
total energy and the energy inside the aforementioned 
shell throughout the simulations. 

While at the early stage there is seemingly linear de- 
crease in energy of the oscillon as shown in Figure [T] this 
rate of energy loss flattens out. Figure [2] shows the energy 
inside the radius r = 5tq over the simulation that reaches 
over 10 s time units. First of all, there is no sign of the 
demise of the oscillon by this point in time. The energy 
is naturally monotonically decreasing function of time, 
but the rate of decrease slows down raising the question 
if the energy approaches asymptotically a constant value, 
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FIG. 2: The energy inside a shell of radius r = 5ro in (j> 4 the- 
ory in two dimensions over the span of the simulation up to 
10 s time units. 
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FIG. 3: The difference of the energy inside a shell of radius 
r — 5ro and the constant Eq = 4.613 in <f> 4 theory in two 
dimensions. The grey straight line is a guide to eye of a power 
law (E - E ) ~ t~ s with the slope 8 = 0.152. 



which wc denote by Eq from now on. We searched for 
a power law of the form (E(t) — Eq) ~ t~ 5 . In practice 
we performed least-square fit with three free parameters, 
the asymptotic energy Eq, the exponent S and an ampli- 
tude in the power law for the data. The best fit for the 
data presented in Figure [2] is shown in Figure [3] where 
we plot the difference of the energy inside the shell and 
the constant Eq = 4.613 on a logarithmic scale. The 
grey line is a guide to eye of a slope —0.152. There is 
very good agreement to this over three orders of magni- 
tude in time providing strong evidence for the power law 
decay of the radiative component governed by a small 
exponent S = 0.1520 ± 0.0005 (we quote hereafter uncer- 
tainties in quantities based on 68% confidence intervals 
derived from standard regression analysis applied to the 
subset of the data used for the fit). Moreover, there exist 
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FIG. 4: The amplitude as a function of time in 4> 4 theory in 
two dimensions. The inset shows the data points in the fit 
A — Ao for Aq = 0.33 for time greater than 10 . The grey 
straight line is the best fit power law (A — Ao) ~ t~ B with the 
slope g — 0.13. 
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FIG. 5: The oscillation frequency as a function of time in 
(j> 4 theory in two dimensions. The inset shows the data points 
in the fit u* for u* = 0.981 for time greater than 10 5 . The 
grey straight line is the best fit power law (oj* — coo) ~ t~ J 
with the slope 7 = 0.182. 

asymptotic value for the energy Eq of the oscillon with a 
value Eq = 4.613 ± 0.003, thus Eq > 4.6. 

The decrease in energy is a consequence of the decline 
in the amplitude of the oscillations, i.e. the maximum 
excursion the field makes at the center of the oscillon, 
which wc take to be towards the local maximum of the 
potential, initially set by C = 1. The amplitude A as a 
function of time is shown in Figure |4] Also this quantity 
is fitted relatively well by a power law, (A(i) — Aq) ~ t~ e . 
The best-fit values for the asymptotic value of the ampli- 
tude is Aq = 0.33 ± 0.06 and for the exponent governing 
the power law g = 0.13 ± 0.02. The comparison between 
the best-fit and the data is shown in the inset in Figure |U 
This suggests that the oscillations would still undergo a 
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FIG. 6: The power spectrum of the field at the center of the 
oscillon, r = 0, in <f> 4 theory in two dimensions. The grey 
straight line is a guide to eye of the exponential exp(— bui) 
with the slope b — 2.2. The inset shows the spectrum around 
the oscillation frequency ujq ~ 0.96. 



substantial decrease in size before reaching the asymp- 
totic value Aq. It should be noted that as the amplitude 
is not an integrated quantity like the energy of the os- 
cillon, the uncertainty in the fit is prominently larger, 
and there remains scatter as can be seen in the inset in 
Figure |U 

With decreasing amplitude the frequency of the oscil- 
lon increases. We measure the frequency through the 
oscillation period, determined from the field's three con- 
secutive crossings of the minimum of the potential at the 
center of the oscillon. Based on the time step used we ex- 
pect the accuracy in determining the frequency from the 
data to be around 0.001. Figure [5] shows the frequency as 
a function of time over the simulation. At the end of the 
simulation the frequency is ui = 0.964, still considerably 
below the threshold for radiation, u/m = 1. 

Previous numerical studies of quasi- or pseudo- 
breathers have established that there exists a critical fre- 
quency lu* such that w* < 1 and above which the os- 
cillon disintegrates when the dimensionality in c/> 4 the- 
ory is more than two (42l . |43| . Furthermore, it has been 
suggested that the approach of the frequency towards its 
critical value is governed by a power law as (lo* — uj(t)) ~ 
i -7 in [43| . Thus we searched again a power law by per- 
forming the least-square fit to the data. The best fit 
yields the critical frequency u* = 0.981 ± 0.003 and the 
exponent governing the power law 7 = 0.18 ± 0.02. The 
data points for uj* — u)(t) are depicted in the inset of 
Figure [SJ There remains a considerable scatter from the 
central value, the grey straight line that demonstrates 
the best- fit power law i -7 , but there is evidence for a 
critical frequency uj* < 1. 

In addition to the oscillation frequency, the whole 
power spectrum of an oscillon is of interest and insight 
to the properties of oscillons can be gained by a study 
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FIG. 7: The Fourier transform of the Jacobi function 
dn(t, 1/2). The height of the first peak is normalised to 10 4 
and its location to unity. 
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FIG. 8: The power spectrum of the field far away from the 
oscillon core at r = 80, in c/> 4 theory in two dimensions. The 
inset shows the spectrum around the location of the most 
prominent peak at u) ~ 1.93. 



in frequency space. We denote the oscillation frequency 
now by ujq to distinguish it from other frequencies present 
and secondly, to highlight that it can be considered a 
constant for any practical purposes in relatively short 
time intervals considered here. A separable ansatz con- 
sisting of a sum of multiplies of this frequency each 
with a separate spatial dependence has seen to converge 
quickly gl IH E| 

^( r ' l ) = ^2 /»( r ) c °s(«W t) . (6) 

71=0 

In an earlier work [4^] we used a method inspired by 
one-particle spectral function at zero momentum in clas- 
sical approximation (for spectral function see [fjfjl l5ll]). 
This technique, employed to trace oscillons moving with 
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varying speed, could be utilised here as well. As the spec- 
tral function involves a volume average it produces gen- 
erally quite a clean signal and reduces background noise. 
However, we wish to study the frequencies present at dif- 
ferent locations on the grid and thus choose to obtain the 
straightforward Fourier transform of the field <p(r) on a 
fixed lattice site over a time interval. 

Figured] shows the Fourier transform of the field at the 
center of the oscillon, (j)(r = 0) over the interval of length 
5000 in time units starting at time t = 1.024 • 10 8 . The 
oscillation frequency wo (thus at this time ujq = 0.964) 
is naturally the most pronounced mode, its multiples ap- 
pearing in the spectra with starkly smaller amplitudes. 
This suppression is clearly exponential, illustrated by the 
grey straight line superimposed above the peaks in Fig- 
ure [SJ One should note here that when performing the 
discrete Fourier transform in a limited time interval it is 
likely that we do not capture the height of the very nar- 
row peaks very precisely. The exponential fit to the data 
with the slope —2.2 should be thus considered rather in- 
dicative than quantitative. The results are qualitatively 
very similar to those obtained in [4~D ] where periodic 
boundary conditions were utilised in a two-dimensional 
lattice. However, the signal here is very clean; there is 
hardly any noise present, but the curve shown is a smooth 
line. We point out that fast Fourier transform of Jacobi 
functions in an interval yields a very similar signal, in- 
cluding the exponential suppression of the amplitudes at 
which the higher multiples of the basic frequency appear 
as well as similar patterns between the peaks. An exam- 
ple is given in Figure [7] where Jacobi function dn (w, 1/2) 
is shown. The resemblance is not co-incidental - Jacobi 
functions are solutions to a type of elliptic non-linear dif- 
ferential equations that result also by omitting the spa- 
tial dependence in the equation of motion @ in </> 4 the- 
ory. Elliptic partial differential equations have also been 
shown to describe quasi-breathers in <f) 4 theory (Tlj . 

The use of absorbing boundary on the grid offers the 
oppourtunity to examine the outgoing radiation far away 
from the oscillon and circumvent the influence of the fre- 
quency directly related to oscillations at u) = ujq. Figurc[5] 
shows the Fourier transform of the field in the same time 
interval as in Figure |H but at location r = 80. This 
choice guarantees that the observation point is a consid- 
erable distance away from the oscillon core, but not at 
the edge of the lattice either where the equation of mo- 
tion j2]) is modified to create the absorbing boundary. 
There is signal at radiation frequency lo — 1 and just 
above it, but the most prominent peak is at the frequency 
lu ~ 1.93, thus coinciding with the second frequency 2wo 
in the expansion ([6]) . There is a peak in the spectra also 
at 3wo and just a visible structure at 4wo, but the signal 
at to = 2wo is the most dominant, the height of the peak 
being over thirty times larger than that at other frequen- 
cies, though still almost six orders of magnitude below 
the main peak at the oscillon core. This observation is 
evidence for that it is the first mode in © for which the 
frequency is above the radiation threshold (here n = 2) 
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FIG. 9: The energy inside a shell of radius r = 5r in (j> 4 the- 
ory in three dimensions. 
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FIG. 10: The difference of the energy inside a shell of radius 
r = 5ro and the constant 41.56 in (f> 4 theory in three dimen- 
sions. The grey straight line is a guide to eye of a power law 
given by (E — Eq) ~ t~ s with the slope 5 = 0.65. 



that is primarily responsible for carrying away the energy 
from the oscillon. 



B. Properties in 4> 4 theory in three dimensions 

In three dimensions we start with a slightly broader 
width but considerably more energetic initial profile than 
in two dimensions as we set C = — 1 and thus the ampli- 
tude of the excursion away from the minimum is towards 
the steeper side of the potential. The energy inside the 
radius r = 5ro throughout the simulation is shown in Fig- 
ure [HI There is a very abrupt drop in energy at the very 
beginning, but then the decrease slows down and energy 
stabilises above the value E sa 40 for a considerable pe- 
riod of time, until the oscillon disintegrates around time 
t = 7200 which appears in the data as the second stark 
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FIG. 11: The amplitude as a function of time in 4 theory in 
three dimensions. The grey line is a guide to eye of a power 
law given by (A-A ) ~ t~ e with A = -1.374 and g = 0.066. 
The inset shows the residuals in the afore mentioned fit. 
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FIG. 13: The power spectrum of the field at the center of the 
oscillon, r = 0, in <f> 4 theory in three dimensions. The grey 
straight line is a guide to eye of the exponential exp( — buS) 
with the slope b = 1.8. 
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FIG. 12: The oscillation frequency as a function of time in 
4> A theory in three dimensions. The inset shows the data 
points in the fit uj* — u> for co* = 1.030 for time greater than 
100. The grey straight line is a guide to eye of a power law 
(lj* — lj) ~ t ' with the slope 7 = 0.2. 

drop in energy. 

Even though the oscillon clearly has a finite lifetime 
here, the steady decrease in energy is relatively well cap- 
tured by a power law. This is demonstrated in Fig- 
ure [TU] where the difference between the energy inside 
the shell and a constant Eq = 41.56 is shown on a 
logarithmic scale. The best-fit value for the exponent 
is S = 0.65 ± 0.01, thus the strength of the emitted ra- 
diation decays far faster in time than in two dimensions. 
This is naturally only of secondary importance because 
the asymptotic value for the energy E = 41.56 ± 0.05 is 
well below the value when the oscillon destabilises, which 
is approximately at E = 42.7. 

The evolution of the amplitude A of the oscillations 
is shown in Figure 1111 On top of the general downward 



FIG. 14: The power spectrum of the field far away from the 
oscillon core at r = 80, in cf> 4 theory in three dimensions. 



trend that follows closely the time evolution of the en- 
ergy, the most notable feature is the growing fluctuation 
in this quantity in the course of time. This amplified os- 
cillation in the amplitude around the central value with 
a decreasing frequency makes a simple power law fit to 
the data challenging. The best-fit to the data in the in- 
terval 500 < time < 5000 yields A = -1.4 ± 0.3 and 
g = 0.066 ± 0.007, but it should be emphasized that the 
values in the fit are very sensitive to the selection of the 
time interval. The residuals of the fit are shown on a lin- 
ear scale in the inset in Figure [TT1 This highlights further 
the growing fluctuations, the beat, in the amplitude; to 
capture the time evolution of the amplitude entirely and 
adequately a model with extra parameters to account for 
these oscillations would be required. The negative value 
of the asymptotic amplitude Aq indicates strongly the 
finite lifetime of the oscillon. 
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The oscillation frequency as a function of time is shown 
in FigureQ21 The frequency approaches u ~ 0.97 towards 
the end of the simulation before the oscillon decays. The 
fluctuations in the amplitude are naturally imprinted into 
the frequency. However, these do not prevent reaching a 
reasonable power law fit to the data over the whole span 
of the oscillon life time. In the best-fit the asymptotic 
value of the frequency is uj* = 1.030 ± 0.003 and the 
exponent governing the power law 7 = 0.20 ± 0.01. The 
power law is demonstrated in the inset in Figure [T2"1 Like 
in the case of the amplitude, the asymptotic value of the 
frequency uj* remains physically meaningless as w* > 1; 
the oscillon has a finite lifetime and decays even before 
the frequency reaches the threshold for radiation uj = 1. 
For more extensive studies of oscillons in quartic theory 
in three dimensions, in particular the effect of the initial 
radius and amplitude, see (45j . 

While the decay of the oscillon in three dimensions re- 
stricts the length of the time interval we can observe the 
evolution, this fate merely enhances the importance to 
examine the frequency spectrum of the oscillon before 
its demise. Due to the relatively fast energy loss the os- 
cillation frequency shifts steadily and to overcome this 
effect we perform the Fourier transform over the inter- 
val of length 200 in time units (the choice of the shorter 
interval naturally compromises the resolution we can ex- 
pect on the frequency axis compared with the analysis in 
two dimensions) . We choose to start the interval at time 
t = 4000, thus this interval is in the middle of the phase 
of relative subdued radiation appearing as the plateau in 
energy in Figure [9j 

The power spectrum at the center of the oscillon, r = 0, 
is shown in Figure 1131 It is qualitatively similar to the 
spectrum in two dimensions (Figure |6|) with peaks lo- 
cated at multiples of the oscillation frequency, which is 
here located approximately at ujq ~ 0.955. There are 
even similar patterns between the peaks as those seen 
in two dimensions. In addition, the suppression of the 
amplitudes which the peaks appear is still exponential as 
well, the best fit value of the slope being —1.8, very com- 
parable with the value obtained in the two-dimensional 
theory. Unsurprisingly, the form of the power spectrum is 
dictated by the potential rather than the dimensionality 
of the theory. 

The major difference between two and three dimen- 
sions is the far larger emitted radiation taking place in 
the latter resulting to a greatly larger width of the peaks 
in frequency space. This has been considered analytically 
in [HI, |45| in the case of the oscillation frequency. 

Figure Q3] shows the Fourier transform of the field c\> 
obtained in the same time interval as in Figure 1131 but 
at the location r = 80. The strength of radiation shows 
simply in the far larger absolute values of the amplitudes 
compared with the situation in two dimensions, see Fig- 
ure [8] There are three distinctive peaks present in the 
spectrum. There is a substantial structure rising imme- 
diately above the threshold for radiation at uj = 1. The 
origin of it can be traced to the large width of the peak 



of the basic oscillation frequency in Figure HH the part of 
this peak for which uj > 1 overflows in the form of free ra- 
diation (for a quantitative study, see [4f| ) . The radiation 
signal at uj ~ 1 is far stronger than in two dimensions. 
However, also in three dimensions the dominant peak is 
at the frequency corresponding the double of the oscil- 
lation frequency, 2wo- The height of this peak is several 
factors larger than the one at uj = 1. Moreover, as waves 
moving at frequency uj = 2ujq carry far larger amount of 
energy than those barely in the free, propagating domain 
at uj > 1, we conclude that the dominant channel for the 
energy loss in three dimensions is also through the fre- 
quency band at uj ~ 2cjo- Finally, there is inferior, but 
still a visible peak in the spectra at uj = 3ujq. 

IV. SINE-GORDON MODEL 

Sine-Gordon model is of considerable interest from the 
point of view of oscillating solutions both on theoretical 
and phcnomenological grounds: in one-dimension it has a 
stable breather solution and it is also the potential of ax- 
ion field (formation of oscillating energ y c oncentrations 
in axion field were studied in |52j). In [53j a fine-tuned 
initial ansatz was reported to create almost dissipation- 
less breather in two dimensions. Here we start with the 
same Gaussian ansatz as in (f> A theory given by ([5]) with 
r ~ 2.9 and C = 1. We scale the potential of sine- 
Gordon model 

y(0) = ^(l-cos(a0)) (7) 

to 

7(0) = ^(l + cos(7r0)) > (8) 

so that locations of the minima are the same as in (f) A 
theory, but here the mass m 2 = 1. 

A. Properties in two dimensions 

Figure [TBI shows the energy inside the shell with radius 
r = 5ro in the simulation up to 5 • 10 7 time units. This 
is qualitatively very similar to the case in quartic the- 
ory. We sought again for a power law using least-square 
fit. The result is presented in Figure [1751 where the con- 
stant Eq = 4.769 is subtracted from the energy inside 
the shell. The resulting curve has the slope —0.119 on 
the logarithmic scale, demonstrated by the grey straight 
line in the figure. Wc lack some of the dynamical range 
reached in the simulation in <fi A theory, but there is evi- 
dence for a power law roughly over three orders of mag- 
nitude in time with an exponent 8 = 0.119 ± 0.003. The 
rate of decrease is marginally smaller than what we ob- 
served in <p A theory, and the asymptotic energy of the 
oscillon, Eq = 4.77±0.01, is now slightly larger, but gen- 
erally oscillons in both models are comparatively similar 
lumps of energy 
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FIG. 15: The energy inside a shell of radius r = 5ro in sine- 
Gordon model in two dimensions over the span of the simu- 
lation up to 5 • 10 7 time units. 
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FIG. 16: The difference of the energy inside a shell of radius 
r — 5ro and the constant Eq — 4.769 in sine-Gordon model 
in two dimensions. The grey straight line is a guide to eye of 
a power law (E — Eg) ~ t~ s with the slope 8 = 0.119. 



The amplitude of the oscillations behaves also very 
similar way compared with the quartic theory in two di- 
mensions. Figure [T7] shows the amplitude over the course 
of time. The inset demonstrates the best-fit power law 
where A = 0.37 ± 0.10 and g = 0.13 ± 0.05. The de- 
rived errors in the fit are quite large, it can be seen that 
there is substantial scatter around the central value for 
time < 10 6 , the oscillon in sine-Gordon model has a larger 
asymptotic amplitude Aq than in 4 theory. 

The evolution of the frequency as a function of time is 
shown in Figure [TH] Like in the case of energy, the rate of 
change is less than in quartic theory and the oscillation 
frequency reaches value uj = 0.961 at the end of the simu- 
lation at time point t = 5.12- 10 7 . The search for a critical 
frequency and an associated power law here yields quite 
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FIG. 17: The amplitude as a function of time in sine-Gordon 
model in two dimensions. The inset shows the data points in 
the fit A - A for A = 0.37 for time greater than 10 5 . The 
grey straight line is the best fit power law (A — Ao) ~ t~ e 
with the slope q = 0.13. 
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FIG. 18: The oscillation frequency as a function of time in 
sine-Gordon model in two dimensions. The inset shows the 
data points in the fit uj* — u for uj* — 0.991 for time greater 
than 10 s . The grey straight line is the best fit power law 
(uj* — Wo) ~ t~ n with the slope 7 = 0.122. 



a stringent results. There is a remarkably good fit to the 
data with uj* = 0.991 ±0.002 and 7 = 0.120 ±0.004. The 
data points for uj* — u){t) are shown in the inset of Fig- 
ure [TH] together with the grey straight line that demon- 
strates the power law i~ 7 . This suggests now that there 
exists a critical frequency w* < 1 in two-dimensional sine- 
Gordon model that would correspond an oscillon with a 
minimum energy 

The power spectrum of the oscillon in sine-Gordon 
model is shown in Figure [T9l obtained by a Fourier trans- 
form in an interval of length 5000 in time units starting 
at time t = 5.12 ■ 10 7 . There are no even multiples of the 
oscillation frequency ujq present in the spectra. This is 
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FIG. 19: The power spectrum of the field at the center of the 
oscillon, r — 0, in sine-Gordon model in two dimensions. The 
grey straight line is a guide to eye of the exponential exp(— b oj) 
with the slope b — 2.1. The inset shows the spectrum around 
the oscillation frequency too — 0.96. 
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FIG. 20: The power spectrum of the field far away from the 
oscillon core at r = 80, in sine-Gordon model in two dimen- 
sions. The inset shows the spectrum around the location of 
the most prominent peak at u — 2.88. 



simply due to the symmetry of the potential ([7|) in sine- 
Gordon model with respect to its minima - there cannot 
be non-zero terms f n {f) in @ for an even n. The remain- 
ing peaks are exponentially suppressed as a function of 
frequency oj. This suppression is even quantitatively very 
similar to the one in quartic theory. The small difference 
in the slope we measure is well within the limits of un- 
certainty we can expect in evaluating the height of the 
narrow peaks. 

Figure [2UJ shows the Fourier transform of the field in 
the same time interval as Figure [19j but at location r = 
80. Now the most dominant peak is at w = 2.88 = 3cj , 
the first radiative mode as the even multiples of ojq, in- 
cluding 2ojq, do not exist. There is still a visible pinnacle 
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FIG. 21: Potential V((j>) in (|10|l around the minimum for 
values p = 0.4, 0.3, 0.2 (from top to bottom). They grey line 
shows the quadratic potential <f> /2. 
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FIG. 22: The second derivative of the potential V{4>) in (fT0)l 
for p = 0.4,0.3,0.2 (from broader to narrower spikes). The 
grey line shows the constant m 2 = 1 for a quadratic potential. 



at 5 ojq and a substantial peak just above the radiation 
frequency oj = 1, but exactly like in cf> 4 theory that one 
is more than order of magnitude suppressed compared 
with the prominent one at the frequency 3ojq- Also in 
sine-Gordon model it is the first non-zero multiple of the 
basic frequency ojq above the threshold for radiation that 
leaks the energy from the oscillon. 

The Gaussian ansatz ([5]) does not create long-lived os- 
cillons in sine-Gordon model in three dimensions. For an 
investigation with fine-tuned initial conditions see (53j . 

V. A CONVEX POTENTIAL 

Majority of the studies on oscillons has been carried 
out in the context of (j) 4 theory (see e.g. 0,0,111,113) or 
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in sine-Gordon model. Both potentials ([3]) and (J7J have 
degenerate vacua and inflections points. In [54[ the study 
dealt potentials that are nearly quadratic. In particular 
a quadratic potential with a negative logarithmic cor- 
rection was reported to support long-lived lumps (called 
I-balls by the authors). Very recently there has been in- 
terest in convex potentials and a study [3lT | considered 
the formation of oscillons in the following potential 



V{4>) 



m 2 M 2 



1 



M 2 J 







1 



(9) 



which is motivated by a number of supergravity and su- 
perstring scenarios. 

It was the work in [54| that initially lead us to study 
the effects of the form of potential on the existence of 
long-lived oscillating solutions. We have chosen to study 
the following convex function 



V{4>) 



l 



(10) 



which reduces to quadratic when exponent p = 0. Fig- 
ure (2JJ shows the potential (fTU)) close to its minimum for 
values p = 0.4, 0.3 and 0.2 together with the quadratic 
potential (j) 2 /2 indicated by the grey line. The potential 
is steeper than quadratic one in the vicinity of the mini- 
mum where oscillations of the field take place, \cf>\ < 1. 

The potential (fTU)) sets far greater challenge for the 
absorbing boundary conditions than <p 4 theory or sine- 
Gordon model. This becomes readily apparent in Fig- 
ure [22] where the second derivative of the potential is 
shown for the same values of parameter p. Already small 
deviations from the vacuum (f> = create large varia- 
tions in V" , the greater the smaller parameter p is. This 
is a drawback from the point of view of the absorbing 
boundary conditions used because they are based on the 
assumption that the potential can be linearised in the 
equation of motion. Once the fluctuations create large 
corrections, the absorption can be expected to be only 
partial or, in the worst case, the update on the boundary 
even pumps energy into the lattice. We restrict consider- 
ations from now on to values 1/2 > p > 0.4 which yield 
less stark variations and for which the radiation is clearly 
adequately removed from the lattice. Potential ([TO]) is 
also more time consuming to evaluate numerically and 
consequently the simulations do not reach such a long fi- 
nal time as the one presented in <fi A theory, but typically 
order of 10 6 time units. 



A. Properties in two dimensions 

Also here the oscillon is initialised by a Gaussian de- 
viation from the vacuum 



4>{r) = C-eM-r 2 /rl), 



(11) 



with the same width as before tq ~ 2.9, but due to the 
steepness of the potential (fTTj]) we set the amplitude of 
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FIG. 23: The difference of the energy inside a shell of radius 
r = 5ro and the constant Eq — 3.80343 in the convex poten- 



tial gOll for p 
line is a 
slope S = 0.97. 



0.45 in two dimensions. The grey straight 
line is a guide to eye of a power law (E — Eq) ~ t~ s with the 



the displacement to be smaller, C = 0.6. The oscillons 
formed are well described by a Gaussian shape though 
with a slightly larger width and a longer tail. After a 
short initial phase of radiation the oscillon settles into a 
state characterised by a very constant energy. We present 
the data from simulation when the exponent was set to 
p = 0.45. The energy inside the radius r = 5ro for this 
choice remains above E = 3.8. We carried out the least- 
square fit for the energy. Now far more fine-tuning is re- 
quired for seeking the asymptotic energy Eq. Figure [23] 
shows the constant Eq = 3.80343 subtracted from the 
energy inside the shell. The radiative component decays 
now almost inversely proportional to time, the best fit 
exponent being S = 0.97 and the corresponding power 
law demonstrated by the grey line in the figure. There is 
agreement with this power law over two orders of mag- 
nitude in time during which the difference (E(t) — Eq) 
drops below 10~ 4 . The best-fit values for the asymp- 
totic energy and the exponent governing the power law 
are Eq = 3.80343(4) and 5 = 0.974 ± 0.003, respectively. 
The diminutive uncertainty in the value of energy reflects 
the fine-tuning needed in the fit, and in turn the stability 
of the value the energy reaches in the simulation. 

Nearly constant energy implies also very stable oscil- 
lation frequency. After initial phase the oscillation fre- 
quency settles to ujq = 0.879 and we do not observe any 
increase in this, but within the accuracy to determine 
frequencies it remains the same until the end of the sim- 
ulation beyond time t = 4 • 10 6 . It is noteworthy that 
this value is considerably below the threshold for radia- 
tion. The amplitude of the oscillations settles to a value 
A ~ 0.49. 

We show the power spectrum of the oscillon in the po- 
tential (|10p for p = 0.45 obtained in the previously used 
time interval starting at time t = 4 • 10 6 in Figure [2U 
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FIG. 24: The power spectrum of the field at the center of the 
oscillon, r = 0, in the convex potential (|10[) for p = 0.45 in 
two dimensions. The grey straight line is a guide to eye of the 
power law with the slope n = 4.6. The inset shows the 
spectrum around the oscillation frequency luq — 0.88. 



10' 




10 ' ' ' 1 ' 

10 2 10 3 10 4 10 5 10 6 10 7 

Time 

FIG. 26: The difference of the total energy and the con- 
stant Eo = 13.908 in the convex potential (|10[) for p = 0.45 
in three dimensions. The grey straight line is a guide to eye 
of a power law (E — Eo) ~ t~ s with the slope 5 = 1.05. 
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FIG. 25: The power spectrum of the field far away from the 
oscillon core at r — 80, in the convex potential (|10|) for p = 
0.45 in two dimensions. The inset shows the spectrum around 
the most prominent peak at frequency uj ~ 1.07. 



Because the potential (fTCTj) is symmetric around its mini- 
mum, there are only peaks at odd multiples of the oscil- 
lation frequency as in sine-Gordon model. But it is sur- 
prising that the suppression is governed by a power law, 
we plot also the frequency axis on a logarithmic scale in 
Figure [24j On the basis of the stability and weak radia- 
tion from the oscillon a strong exponential suppression of 
radiative modes would have been likely. The power law 
decrease as a function of the frequency is considerably 
quick though, the amplitudes are suppressed roughly as 
lu~ 5 though the exponent should be considered only in- 
dicative. 

Figure [51] shows the frequencies present far away from 



the oscillon core at r = 80. There is a peak at fre- 
quency lu = 2.64, thus at 3wo, but this one is not the 
dominant. Instead the amplitude of the peak just above 
the radiation frequency at w = 1.07 is over seventy times 
higher. The dominant radiative mode is not related to 
multiples of the oscillation frequency luq. 

B. Properties in three dimensions 

We use the same Gaussian initial ansatz (|11[) with the 
same width and amplitude as in two dimensions. This 
initial profile creates now an oscillating lump also in three 
dimensions though the oscillon has a larger width and 
a longer tail than the Gaussian distribution. The data 
shown is for the same value p — 0.45 in the potential (|10|) . 
After initial stronger radiative phase, there remains en- 
ergy above the value E = 13.9 within the radius r = hr$. 

For the least-square fit we use now the data on the total 
energy in the lattice because the oscillon is more spatially 
spread and due to the precision required in determining 
the value of the constant Eq below 10" 3 the fluctuations 
in the energy inside the shell become apparent (they are 
to a lesser extent already visible in Figure |2"3")) . Figure l2"6l 
shows the constant Eq = 13.908 subtracted from the to- 
tal energy. The decaying component decreases fast in 
time, the best-fit for the exponent governing the power 
law is 6 = 1.09 ± 0.01 while the asymptotic value of the 
energy is E = 13.9081 ± 0.0002. There is fairly large 
uncertainty in the value of the exponent, but there is 
good evidence that the radiative part decays faster than 
inversely proportional to time. The amplitude related to 
this energy is A ~ 0.34. 

With the aforementioned energy the oscillation fre- 
quency settles to lo ~ 0.939. Again we do not observe 
increase in the frequency within the numerical accuracy 
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FIG. 27: The power spectrum of the field at the center of the 
oscillon, r = 0, in the potential (|10|l for p = 0.45 in three 
dimensions. The grey straight line is a guide to eye of the 
power law uj~ k with the slope k = 4.85. The inset shows the 
spectrum around the oscillation frequency uio ~ 0.94. 
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FIG. 28: The power spectrum of the field far away from the 
oscillon core at r = 80, in the potential (| 10[l for p = 0.45 in 
three dimensions. The inset shows the spectrum around the 
most prominent peak at frequency uj ~ 1.04. 



tion frequency. The only peak of any note is located at 
frequency u> = 1.04, thus barely above the threshold for 
propagating radiation. 

To summarize our findings of the oscillons in the con- 
vex potential (|10[) . the change from two to three dimen- 
sions seem to have surprisingly small effect. In both di- 
mensions there remains a very stable lump of energy os- 
cillating at a constant frequency which emits only tiny 
amount of radiation. While we do not have an expla- 
nation for their longevity, the investigation in frequency 
space gives one answer. As the oscillon can primarily 
only excite modes that are just above the threshold for 
radiation, i.e. (cj — m) -C 1, obviously such modes can 
carry only lesser amount of energy away from the core of 
the lump. Consequently we can expect far reduced emit- 
ted radiation compared with c/> 4 theory and sine-Gordon 
model as we indeed observe (it cannot be excluded that 
some of the signal we see is from relic radiation on the 
lattice that was not absorbed by the boundary, but re- 
flected back to the grid). 

We have considered here only very limited range of 
values of the parameter p in the potential (flTj|) , but the 
potential supports oscillons in a wider variet y o f the val- 
ues of the exponent p governing its steepness [5 51 ] . With a 
different choice for lattice spacing and time step, oscillons 
have survived up to 10 7 time units [5oT ]. This is hardly 
a surprising result taking into account their minuscule 
radiation losses. 

Finally, it should be noted that non-radiating oscil- 
lons has been found in one-dimensional signum-Gordon 
model [56[ (more recently also a new class of swaying os- 
cillons have been discovered in the same model [57j). It 
is defined by the potential 

V{ct>)=g\ct>\. (12) 

The convex potential (fTTJ)) reduces to linear for p = 1/2 
in the limit of large field, \<p\ 1. This is not the realm 
where oscillations of the field take place, but it cannot 
be excluded right away that there were no connection 
between oscillons we have reported here and those de- 
scribed in (56|. At very least, both are examples of os- 
cillons appearing in theories where the potential is not a 
smooth one. 



after the initial phase. Figure |2"T1 shows the power spec- 
trum of the three-dimensional oscillon obtained in the 
previously used time interval starting at time t = 4 • 10 6 . 
The result is very similar to Figure [M] with peaks at odd 
multiples of the basic frequency loq. The suppression of 
the higher amplitudes is again governed by a power law 
demonstrated by the grey straight line in the figure. The 
measured slope is slightly steeper than in two dimen- 
sions, but this difference is hardly significant taking into 
account the accuracy of the data we can expect. 

Figure |2"81 shows the frequencies at the location r = 80. 
There is only barely visible structure above the back- 
ground at u> ~ 2.82, three times the value of the oscilla- 



VI. CONCLUSIONS 

We have studied oscillons in <f> A theory, sine-Gordon 
model and in convex potentials, performing (1+1)- 
dimensional simulations assuming radial symmetry. An 
important ingredient was the development of absorbing 
boundary conditions for a massive, real scalar field. Gen- 
erally absorbing boundary conditions work well and ra- 
diation is removed from the lattice. 

We showed that the time evolution of the energy and 
the frequency of oscillons is well modelled by dividing it 
to a constant and a decaying radiative part. There is 
strong evidence that the decay of the radiative part is 
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model 


dim 


E 


8 


quartic theory 


2 


4.613 ± 0.003 


0.1520 ± 0.0005 


quartic theory 


3 


41.56 ±0.05 


0.65 ±0.01 


sine-Gordon 


2 


4.77 ±0.01 


0.119 ±0.003 


convex potential 


2 


3.803430 ± 4 ■ 10~ b 


0.974 ± 0.003 


convex potential 


3 


13.9081 ± 0.0002 


1.09 ±0.01 



TABLE II: The values of the asymptotic energy Eq and the 
exponent S in (E(t) — Eo) ~ t~ s for the oscillons in the var- 
ious models studied with the initial conditions as mentioned 
in the text. 



model 


dimension 


^0 


Q 


quartic theory 


2 


0.33 ± 0.06 


0.13 ±0.02 


quartic theory 


3 


-1.4 ±0.3 


0.066 ± 0.007 


sine- Gordon 


2 


0.37 ±0.10 


0.13 ±0.05 



TABLE III: The values of the asymptotic amplitude Ao and 
the exponent g in (A(t) — Ao) ~ t~ e for the oscillons in the 
various models studied with the initial conditions as men- 
tioned in the text. 

governed by a power law in all models considered. We 
further determined the values of the constant components 
as well as values of the exponents in the power law decay, 
which we summarize in Table |TTJ The power law decay 
has been previously reported in 4 theory [H, [45J and 
for flat-top oscillons in [20j . Note that for the oscillons 
in the convex potentials the frequency remains constant, 
which is a qualitatively new behaviour for an oscillon. 

We also studied the evolution of the amplitude and the 
oscillation frequency over time, our results are summa- 
rized in Tables IIIII and IIVI The oscillon amplitude also 
exhibits a power law decay, which in three dimensions 
is strongly modulated by a beat frequency, while in the 
case of the convex potential it is constant. 

Furthermore, we examined the power spectrum of field 
at the centre of oscillon and near the boundary. While 
in </> 4 theory and in sine-Gordon model there is an expo- 
nential attenuation of higher multiples of the oscillation 
frequency, the oscillons in the convex potential show only 
a power law suppression. 

Understanding the manner oscillons radiate and shed 
their energy away is crucial in gaining insight for their 
persistence and longevity. We studied the outgoing ra- 
diation in frequency space. The emitted radiation in </> 4 
theory and in sine-Gordon model is strongly peaked at 
the first non-zero multiple of the basic frequency, while 
in the convex potential the oscillons radiate at a low fre- 
quency just above the threshold to release radiation. Our 
results in quartic theory, in particular in two dimensions, 



model 




7 


quartic theory 


0.981 ± 0.003 


0.18 ±0.02 


sine-Gordon model 


0.991 ± 0.002 


0.120 ±0.004 



TABLE IV: The values of the critical frequency to* and the 
exponent 7 in (uj* — oj(t)) ~ i -7 in two dimensions for the 
oscillons with the initial conditions as mentioned in the text. 



agree well with the study of quasi breathers in the same 
model |43j . There it was identified that the emissions 
take place at the frequencies nuiQ, with n > 2 with the fre- 
quency 2 wo being the dominant one. The quasi breather 
approach, i.e. ansatz ([6]), does not assign any width for 
the frequencies. This was considered the starting point in 
the work (45j when dealing with the basic oscillation fre- 
quency. Our results indicate that there is a contribution 
from both sources, from higher modes as well as from 
the spread related to the basic oscillation frequency and 
once again highlight the richness and complexity of os- 
cillons. An approach that would account for both higher 
frequencies present in the spectra as well as the width 
of the frequencies might capture the problem of oscillon 
radiation in full. 

What do the results tell about the lifetime of oscil- 
lons? The appearance of the constant Eq now raises the 
question of its physical interpretation. While it would 
be tempting to define it as the energy of a non-radiating 
oscillon, this kind of conclusion cannot be made straight- 
forwardly. It is namely not obvious at all that there 
are oscillon solutions at such energy. While the energy 
decreases the oscillon oscillates faster and the oscilla- 
tion frequency u>q increases. Once close enough to the 
threshold for radiation u> = m oscillon demises when it 
reaches its minimum energy state [42T [43j . This was well 
demonstrated in 4 theory in three dimensions, there the 
asymptotic energy Eq is well below the energy level that 
caused the oscillon to dissolve. Because the radiative 
part shows a slow decay in time in two dimensions, the 
energy oscillons have reached at the end of simulations 
is still well above the asymptotic value Eq in 4 theory 
and sine-Gordon model, more precisely by 13% and 7%, 
respectively (it should be noted that there are some the- 
oretical bounds on the minimum energy of lumps [58| ) . 

However, the case for the convex potentials presented 
here is definitely intriguing. We have shown that there 
the suppression of the radiation losses is quick, roughly 
inversely proportional to time. Thus only very small fur- 
ther decrease in energy is to be expected to occur, orders 
of magnitude less than in 4> 4 theory or in sine-Gordon 
model. Furthermore, the frequency of oscillons is well 
below the threshold for radiation. This all makes it pos- 
sible that these oscillons could go on indefinitely long. 
Another open question remaining is if it is possible to 
tailor potentials in which the emission of radiation from 
oscillon is strongly suppressed also in higher dimensions. 
The convex potential we considered is not phenomeno- 
logically particularly well motivated, but if this sort of 
oscillons can also appear in more realistic theories their 
long life-time would almost certainly make them viable 
to affect physics in the system there. 



Acknowledgments 

P. S. thanks Nick Manton, Fuminobu Takahashi and 
Anders Tranbcrg for discussions and is grateful to Melvin 



15 



Vamghese for comments on regression analysis. The au- 
thors acknowledge Robert Bennett for sharing his results 
with them. P. S. is supported by the Claude Leon Foun- 
dation. 



the boundary. Expanding the square root to second order 
in (|A4|) and replacing the variable w by the correspond- 
ing time derivative one obtains an absorbing boundary 
condition on the edge of the one-dimensional grid 



Appendix A: Absorbing Boundaries 



d dip 
dt dr 



d 2 ip 



1 2 1 dip 



0. 



(A5) 



In this appendix the absorbing boundary conditions 
are presented. As the starting point consider the follow- 
ing linear wave equation 



<P 



vV 



<p- V\ =0 = 0, 



(Al) 



where ip is now a small perturbation around a chosen 
vacuum and V" evaluated at that point. For example in 



4> theory tj> v — ip and V"{<j>)\^ =v = 2Xv 2 



i 2 . The 



treatment here follows closely the one presented in [59| for 
a massless scalar field (for absorbing boundary conditions 
in general see also [6(| ) ■ 



1. Cylindrical Wave 

When one assumes radial symmetry in two space di- 
mensions the wave equation (|A1| reads 



dt 2 dr 2 r dr 



(A2) 



The solutions of (|A2p are Bessel functions, but consider- 
ing an absorbing boundary at far distance (r> 1) the 
approximative solution reads 



<p(t,r) = 



1 



: exp 



i ( \/ ui 2 — m 2 r + tot) 



(A3) 



where w is a dual variable to t. This solution (|A3|) is 
annihilated at the boundary r = rb (to order r -5 / 2 ) by 
the operator in the parenthesis 



d_ 

dr 



1 

2~7 



(A4) 



The equation (|A5j) permits a straightforward time evolu- 
tion using leap frog algorithm. To be more precise, the 
condition (|A5[) yields a new evolution equation for the 
field momentum LI compared to that in the interior of 
the lattice, whereas the update of the field </> stays un- 
altered. The last term in (|A5|) . suppressed by inverse 
of physical distance r, results from the treatment in the 
polar coordinates. It turns out to be fairly insignificant 
with the numerical set-up used in this study, but in a 
smaller grids it increases the absorption. 



2. Spherical Wave in Three Dimensions 

When one assumes spherical symmetry in three dimen- 
sions the solution for the wave equation 



d 2 ip d 2 ip 2 3» j 
dt 2 - dr z r dr 



is given by 



ip (t, r) = — exp i (\J lo 2 — m 2 r + uit ) 



This is annihilated at the boundary r = n, by 



i^/lo 2 - 

dr 
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(A6) 



(A7) 



(A8) 



which leads to following absorbing boundary condition 



d dip d 2 ip 1 



1 dip 



di^ + W-2 m{p+ rm = ^ 



(A9) 



From the point of view of numerical implementation us- 
ing leap frog update we aim to second order equations at 



appropriate when a three-dimensional system is consid- 
ered in a spherically symmetric geometry. 
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